a = [];
for alpha=0.01:0.01:0.99
    z=sqrt(3)/pi*(log(alpha)-log(1-alpha));
    y=(exp(1+z)-1)/(1+z);
    a= [a;y];
end
plot(a,0.01:0.01:0.99)
